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Abstract 

We report in details the observations of structures in coupled map lattice dur- 
ing its chaotic evolution, both in one and two dimension, driven by identical 
noise on each site (by a structure we mean a group of neighboring lattice-sites 
for whom values of dynamical variable follow certain prespecified pattern). It 
is observed that the number of these structures decays with their size fol- 
lowing a power-law, for a given noise-strength. The number of structures 
decays with their lifetime following a stretched exponential. We have seen 
an interesting phenomenon, which we call stochastic coherence, in which the 
average length as well as the average lifetime of these structures exhibit bell- 
shaped maxima for some intermediate range of noise-strength values. Similar 
features have also been observed for 'spatio-temporal structures'. 
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I. INTRODUCTION 



Over years coupled map lattice (CML) has become a very popular model to study a host 
of physical phenomena involving systems with spatial dimension [PH^. In particular, CML 
exhibits a variety of spatio-temporal patterns and structures and has been extensively used, 
with varying degree of success, to model similar features in experimental systems such as 
Rayleigh-Benard convection, Taylor-Couette flow, B Z reaction etc. [|r|MB]H|nj . One of the 
challenging problems is to understand the formation of structures, localized in both space 
and time, in turbulent fluid [|I5|-|TT||. Recently lot of attention has been devoted to the role 
of fluctuations in onset, selection and evolution of such patterns and structures |TB|-f2^]. An 
interesting and counterintuitive observation is that the presence of noise can help sustain 
structures which otherwise would have been absent in an evolving dynamical system. 

In this paper we present the detailed account of a novel phenomenon that we have 



recently observed |27] in the dynamics of structures in a chaotically evolving CML, in one 
and two dimension, driven by identical noise. We call a structure as a group of neighboring 
sites whose variable-values follow certain prespecified spatial pattern. As the lattice evolves 
under the influence of noise, the number of these structures exhibits power-law decay with 
size of the structure for a given noise-strength, with an exponent which is a function of 
noise-strength. It is observed that average size of these structures, plotted against noise- 
strength, shows a bell-shaped curve with a characteristic peak. Similar features have been 
noted for the 'spatio-temporal structures' also (by this we mean the cluster of sites which 
remain in a structure as it evolves in time). Average lifetime of these structures also exhibit 
a maximum within same range of noise values. The behaviors remain essentially the same 
in two dimensional CML as well. We call this new phenomenon stochastic coherence. 



II. THE SYSTEM 



We first consider a one dimensional CML with the dynamics 

x t+1 (i) = (1 - e)F{x t {i)) + E - [F{x t {t - 1)) + F{x t {i + 1)) 



+ Vt 



(1) 



where x t (i), i = 1,2, ■ • • , L, is value of the variable located at site i at time t, r\ t is the 
additive noise, e is the (nearest neighbor) coupling strength, and L is the size of the lattice. 



Evolution on the lattice sites is governed by the nonlinear logistic function 

F(x) = fix(l -x), (2) 

where /i is the nonlinearity parameter. Both open-boundary conditions Xt(0) = Xt(l), Xt(L+ 
1) = Xt(L), and periodic-boundary conditions x t (L + i) = x t (i), have been used for our 
system. For noise rjt we have chosen a uniformly distributed random number bounded 
between — W and +W. We call W the noise-strength parameter. 

We define a structure as a region of space such that the dynamical variables at sites 



within this region follow a prespecified spatial pattern [13|]. For our purpose we choose a 
spatial pattern where the difference in the values of the variables of neighboring sites within 
the structure is less than a predefined small positive number say 5, i.e., 

\x t {i) - x t {i ± 1)| < 5 . 

We call S the structure parameter. We study the distributions of abundance and lifetimes 
of these coherent structures as the CML evolves chaotically. 



III. RESULTS AND DISCUSSION 



Values of the system parameters /i, e and L are so chosen that the dynamics of the system 
remains essentially chaotic. We have taken typical values L = 10000, fi = 4 and e = 0.6. 
Value of structure parameter 5 is chosen as 5 = 0.0001. Coherent structures with length 
< 3 (sites) and lifetime < 2 (timesteps) are disregarded. W is restricted to within [0, 1]. For 
the results presented here open-boundary conditions are used. 80000 transient steps have 
been discarded and the data are integrated over 100000 iterates per initial condition and 4 
different initial conditions. 



A. Distribution of n(l) vs. I 

We have investigated the distribution of number n(l) of structures with their length I, 
for different values of noise-strength parameter W. In Fig. 1 we show the plot (on log-log 
scale) of n(l) against I, for five values of W. The plot exhibits, over a wide range of scales, 
a characteristic power-law of the form 
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n{l) oc r ai , 



(3) 



where ot.\ is power-law exponent. 

If one uses only a uniform-deviate random number as site- variable, probability that a 
site belongs to a structure of length I is pg(l) ~ /(2<5)' _1 (1 — 25) 2 . So the number of such 
structures in a lattice of size L will be 

n s {l)^L{25) l -\l-25)\ (4) 

which shows an exponential decay. Thus correlations in the system are important to get a 
power law form. 



B. Stochastic coherence 

One can see from Fig. 1 that the exponent a.\ of relation (^) depends on noise-strength 
W . Fig. 2 shows the variation of ot\ with W . ot\ is exhibiting a minimum for values of W 
around 0.6. To get a clearer picture, let us define average length I oi a structure as 

r=E*»(0/E»(0- (5) 

We have investigated how noise influences this quantity. In Fig. 3 we plot the variation 
of I with W for values of parameters as in Fig. 1. One sees a quite interesting bell-shaped 
dependence for W around value 0.6. 

This behavior may remind one of the much studied stochastic resonance phe- 
nomenon [p8|-|3"0"|j , a signature of which is the existence of a bell-shaped behavior of the 



temporal response of the system, when plotted against noise-strength. The similarity, how- 
ever, is purely coincidental because there exists one important difference. In stochastic 
resonance noise transfers energy to the system at a given time scale, which usually is the 
characteristic scale of the system. On the other hand, power law feature of Fig. 1 implies 
that our system does not have any preferred length scale; noise induces coherence on all 
length scales. We call this novel phenomenon stochastic coherence. 



4 



C. Stability analysis 



Let us see, qualitatively, how noise modifies the stability matrix of the system, thereby 
influencing its stability properties. We consider stability matrix M of a homogeneous state 
{• • • , x t , x t , x t , • • •} at time t. This state may be imagined as a large structure with 5 = 0. 
One time step later the stability matrix is 



M< 



t+i 



JF' 



(6) 



where J is the familiar tridiagonal matrix 



J = I 



•• l-e e/2 
e/2 l-e e/2 
e/2 l-e 



(7) 



with 1 — e as its diagonal elements and e/2 as offdiagonal elements on either sides, and 



F t ' = F'(:r t )(= —)=^l-2x t ). 
After two more timesteps, the stability matrix takes the form 



(8) 



M t+1 M t+2 M t+3 = J^F' t F' t+1 Fl +2 

= J 3 F//i 2 [(l - 2F t ){\ - 2//F t (l - F t )} + 6/i^(l - 2F t ) 

-2rj t {l + - 6F t + 6F, 2 )} - W - 2^ +1 (l - 2F t ) + 4 TO+1 ] . 

We now average this expression over the noise distribution. The noise has uniform distri- 
bution, with zero mean and delta-correlation. By averageing, the terms containing cross 
terms in 77 (i.e. terms of the type i] t r] t+ i) and odd powers in r\ vanish, leaving the following 
expression with only even powers in r\: 



< M t+1 M t+2 M t+3 > = J 3 (F/) 2 /i 1 - 2/iF,(l - F t ) + 6/1 <^> 



(9) 



where <> denotes averaging over noise-distribution. One more averaging is needed, this 
time over the invariant distribution of the CML. Fig. 4 shows the invariant density of a 



single site of CML, with W = 0.6 and the other system parameters same as mentioned in 
the beginning of this section. The density is highly asymmetric, with larger weightage for 
x(i) > 0.5. Because of this asymmetry, averaging expression (§) over the invariant density 
makes the term 1— 2fiF t (l— F t ) negative. This, added to the positive noise term 6fi<r]f>, 
results in reduction of eigenvalue of the above matrix. This in turn implies a consequent 
reduction of the instability of the state (or the structure). Thus noise plays a crucial role in 
enhancing the abundance of the structures. 

D. Evolution of structures 

Let us now consider the evolutionary aspects of these coherent structurs. We have studied 
the distribution of the number n(r) of structures having lifetime r. In order to obtain r, 
each structure is tracked as the lattice evolves, till the structure degenerates completely. 
In Fig. 5 we show such an evolutionary diagram for a CML with a typical size L = 100, 
evolving chaotically for 200 timesteps. The shaded regions depict the evolution of coherent 
structures. Fig. 6 shows the variation, on log-linear scale, of n(r) with r, for different values 
of noise W. Number n(r) exhibits a decay with r with a stretched exponential form 

n(r) oc exp ( — (const. )t iS ) . (10) 

The exponent (3 is seen to depend on W . Like in case of a of expression @, (3 also exhibits 
a minimum near W = 0.6. 

We define average lifetime f of a structure as 

f = 5>n(r)/I>(T) ) (11) 

and explore its dependence on noise strength W . Fig. 7 shows a plot of f vs. W. It again 
exhibits a maximum for W around 0.6. 

There is another feature of the system which merits an elaborate study. This concerns 
the distribution of 'spatio-temporal structures (STS)'. By such a structure we mean the 
group of sites that remain in a coherent structure as it evolves in time, till it completely 
degenerates. Each shaded region of Fig. 5 constitutes one such STS. In Fig. 8 we plot the 
number n(Si) of an STS vs. its 'size' Si for different noise values W. We again see an 
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interesting power law dependence of the type (|3|) for an appreciable range of sizes, but with 
different exponent. This exponent also depends on W, with a minimum around W = 0.6. 
Let us define average size Si of an STS in a similar fashion as 

S 1 = Y f 3M3 1 )/ y £n{3 1 ). (12) 

We plot Si as function of W in Fig. 9. This plot also exhibits a characteristic bell shape 
within a range around W = 0.6, quite similar to the plot in Fig. 3. 



E. Lyapunov spectrum analysis 

We have studied the lyapunov exponent spectrum (A spectrum) of the CML, and its 
dependence on noise W as well as on the coupling parameter e. We find that a number 
of lyapunov exponents is positive, thereby implying that the system is indeed evolving 
chaotically. In Fig. 10 we plot the variation of maximum exponent A max with noise strength 



W, for coupling parameter e = 0.6. A max exhibits a clear minimum around W = 0.6 ||31|| , the 
same region where all the extrema reported above are located. This observation may create 
the impression that the reduction of lyapunov exponent alone is responsible for stochastic 
coherence. 

To explore this possibility, we have studied the dependence of both A max and average 
length / on coupling strength e. Fig. 11 plots the variation of A max with e for different W. 
It is clear that A max remains fairly constant for 0.2 < e < 0.8 for all W. On the other hand, 
/ shows a monotonic increase with e for all W. In Fig. 12 we show one such plot of I against 
e, for W = 0.6. These two contrary features indicate that lyapunov exponent alone is not 
good enough to characterize the full range of spatio-temporal features of the system. 

It is worth pointing out that stochastic coherence is quite distinctive phenomenon as 



compared to spatio-temporal intermittency [32]. In intermittency the system exhibits al- 



ternately regular and chaotic bursts in time. Our system does not show any such behavior. 
This fact is further corroborated by the fact that the power spectrum of time-series does not 
have any peak for entire range of noise-strength. Thus our system is undergoing essentially 
a spatially intermittent and temporally chaotic evolution. 
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F. Full lattice coherence (FLC) 



One quite interesting phenomenon that we have observed is that for W > 0.4 often the 
entire lattice itself evolves as a single coherent structure. This state should be distinguished 
from the synchronized state where at any given instant t, x t (i) = x t for i — 1, 2, • - • , L (that 
is, 5 = 0). On occasions this full lattice coherence (FLC) survives for quite long durations 
(as long as 200 timesteps or more), but eventually degenerates completely. This implies 
that synchronized state is not a stable attractor for our system. We have confirmed this by 
calculating the largest Lyapunov exponent for the synchronized state which turns out to be 
positive. 

However, occasionally our system did fall into an apparently synchronized state after a 
very large time (larger than 10 7 steps). This happens because of finite accuracy of computa- 
tion which cannot differentiate between unstable and stable synchronized states |33|]. Let us 



investigate this point further. We define a quantity T as the average time required for first 
occurrence of FLC, the averaging taken over different initial conditions. In Fig. 13 we plot, 
on log-log scale, T against the structure parameter 5, for a fixed lattice size L = 10000. The 
plot exhibits a power law of the form 

f L oc (J-t . (13) 

This relation implies that the synchronized state (5 = 0) will not occur in a finite lattice. 
We now show, in Fig. 14, the plot of T against lattice size L, for a fixed 5. This is again a 
power law of the type 

f s oc U . (14) 

We thus see that as L — > oo, even FLC with a nonzero 5 will not be attained. Therefore, 
relations ( |T3"D and (|T4"D tell us that synchronized state for our system is an artifact, resuting 
due to the combined effect of finite lattice size and finite computational accuracy. 

Existence of FLC in a finite lattice can be seen to be a consequence of the power-law 
(0). Probability of a site to belong to a structure of length I, following (^), is ps(l) oc 
l.l~ ai = The probability of a site to belong to a structure of length > L is then 

= Ps(> L) oc El P~ ai « II 'dll 1 '* 1 oc L 2 ~ ai (this holds for a x > 2, which is true for 
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our system as can be seen from Fig. 2). Thus, the probability that for the first time a site 
belongs to a structure of length > L at timestep T is Ps(T) oc (1 — P ) T_1 P - Average of T 
is f s (> L) = Et=iTPs(T) oc (P )-\ i.e., 

f 5 (> L) oc L ai ~ 2 . (15) 

This shows that FLC (with nonzero 5) will occur in a finite lattice. From Fig. 2 we get 
oi\ — 2 « 0.22 for = 0.6. On the other hand, we estimate (from Fig. 14) v of relation 
(|I~4|) to be approximately 0.3. The discrepancy may occur due to the fact that relation ( p~5|) 



holds for an infinite lattice, whereas relation flT4|) is numerically obtained for finite lattices 



IV. TWO DIMENSIONAL CML 

Let us now consider a two dimensional CML. The dynamics of this system takes the 
form 

x t+1 {i,j) = (1 - e)F{x t (i,j)) + ~[F(x t {i - 1, j)) + F{x t {i + 

+F{x t {i,j - 1)) + F(x t (i,j + 1))] + m , (16) 

where x t (i,j), i,j = 1, 2, • • • , L, is the value of the variable at site in the lattice having 
size L x L, and F(x(i,j)) = \i x(i,j)(l — x(i,j)). For structures in this lattice we look for 
those sites with 

\x t (i,j) - x t (i ± l,j)\ < 5 , 

and, 

\x t (i,j) - x t (i,j ±1)| < 5 . 

As the CML evolves, we study the distributions of abundance and lifetimes of these two 
dimensional structures, as we have done for one dimensional system. 

We present some of the observations in the following. Again we have chosen the system 
parameters so that the CML evolves chaotically. Typical values taken are Lx L = 100 x 100, 
/i = 4, and e = 0.6. We have chosen structure parameter 5 = 0.0001. 3000 transient steps 
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have been discarded, and data are obtained over 20000 iterates per initial condition and 4 
initial conditions. We have used open boundary conditions for the results presented here. 

Fig. 15 exhibits the plot (on log-log scale) of distribution of the number n(s) of coherent 
structures against their size s for different W (we use the notation s to denote size of the 
structures in two dimensional CML, unlike I which has been used for one dimensional case. 
This is also to be distinguished from Si, the notation used for sizes of STS). One can see a 
power law of the form 

n(s) oc s~ a2 , (17) 

similar to the relation for a wide range of values of s (for s > 10; there is a slight bend 
in the graphs for s < 10). The exponent a 2 for the power law is slightly less in value than its 
one dimensional counterpart ai, for all W, indicating a less steeper decay. These exponents 
again show a dip around W = 0.6. We define the average size s of the structures as 

S = 5>n(a)/5>(*). (18) 

s exhibits a corresponding peak around W = 0.6. This can be clearly seen in the bell-shaped 
curve of Fig. 16, which is the plot of s against W. This plot is quite alike the one of Fig. 3. 

In Fig. 17 we plot, on log-linear scale, distribution of number n(r) of structures against 
their lifetime r. The plot shows a stretched exponential type behavior of the form (p~0|) . 



similar to Fig. 6. Fig. 18 shows the variation of average lifetime f, as defined in (|Tl~D, with 
noise strength W. We see a peak near W = 0.6, quite similar to the plot of Fig. 7. 

Distribution of number n(5' 2 ) of STS is plotted, on log-log scale, against size S 2 in Fig. 19. 
We see a very similar power law decay as in Fig. 8 for one dimensional CML. As in previous 
instances, we define average size S2 of STS as 

S 2 = J £S 2 n{S 2 )/ J £n(S 2 ). (19) 

Fig. 20 shows the variation of S 2 with noise strength W. The plot shows a bell-shaped peak 
around W = 0.6, similar to the plot of Fig. 9. 

All the other features that are observed in one dimensional system are also seen in two 
dimensional CML. For instance, we have seen full lattice coherence for W > 4, which persist 
for considerable duration before eventually breaking up. 
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We have repeated all the above observations, for both one and two dimensional systems, 
for several other values of the coupling parameter e ranging from 0.1 to 0.9, as well as for 
nonlinearity parameter fi between 3.6 and 4. We find that all the features remain essentially 
the same, indicating the robustness of the phenomenon. We have also observed similar 
results using periodic-boundary conditions for the lattice. 



V. CONCLUSION 



We have investigated in details a new phenomenon associated with structures in a chaot- 
ically evolving CML, in both one and two dimensions, driven by identical noise. We call this 
phenomenon stochastic coherence. We have seen a marked rise in abundance of coherent 
structures of all scales due to noise. Stability matrix of these structures shows that noise 
can reduce their instability, thereby enhancing abundance.. Distribution of the structures 
exhibits a power-law decay with size of the structure, with an exponent having a minimum 
at some intermediate noise-strength. Average size of these structures shows a bell-shaped 
maximum at the same value of noise-strength. This feature is similar to that of the stochas- 
tic resonance phenomenon. However, our system does not have any intrinsic length-scale, 
whereas stochastic resonance is associated with a given time-scale. We have also observed 
similar maxima for the average lifetime of these structures as well as for the average size 
of 'spatio-temporal structures', at the value of noise close to the earlier extrema. Our phe- 
nomena may have their importance in understanding the interesting role that noise plays in 
formation and evolution of structures in spatially extended systems. 
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FIGURE CAPTIONS 



Fig. 1. The figure shows the variation, on log-log scale, of number n{l) of structures with 
length / for a one dimensional lattice with size L = 10000, for different values of noise- 
strength W as indicated. The system parameters are coupling strength parameter 
e = 0.6, structure parameter 5 = 0.0001, and nonlinearity parameter ji — 4. Open- 
boundary conditions are used. 80000 transient steps are discarded. Data are obtained 
for 100000 iterates per initial condition and with 4 initial conditions. 

Fig. 2. Exponent «i of the power law (|j) (slope of the plots in Fig. 1) is plotted against noise 
strength W, with system parameters as in Fig. 1. 

Fig. 3. Variation of average length I of structure with W is plotted, with parameters as stated 
in Fig. 1. Vertical error bars indicate the standard deviations at the data points. 

Fig. 4. Invariant density p(x) of the CML is plotted against variable value x, for the system 
parameters as for Fig. 1 and with noise value W = 0.6. The asymmetry is clearly 
evident, with a much higher weightage for x > 0.5. 

Fig. 5. Evolution of coherent structures is shown by shaded regions, during 200 time steps, 
in a CML with a typical size L = 100, and with noise value W = 0.6, and other 
parameters same as in earlier figures. 

Fig. 6. Distribution of number n(r) of the structures is plotted, on log-linear scale, against 
their lifetime r. The system parameters are as in Fig. 1. 

Fig. 7. Plot shows the variation of average lifetime f of structures with noise-strength W, for 
parameters as in Fig. 1. 

Fig. 8. Distribution of the number n(Si) of spatio-temporal structures (STS) is plotted, on 
log-log scale, against their size Si, for the same parameter values as in Fig. 1. 

Fig. 9. The plot shows the variation of average size Si of STS with the noise-strength W, for 
parameter values as in earlier figures. 
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Fig. 10. Variation of A max is plotted against noise strength W, for e = 0.6. 

Fig. 11. The plot shows variation of A max with coupling strength e, for different values of W, 
other system parameters remaining the same as in Fig. 1. 

Fig. 12. Average length I of the structures is plotted against e, for W = 0.6. 

Fig. 13. Variation of T is plotted, on log-log scale, against S, for given lattice size L = 10000, 
other parameters being as in earlier figures. 

Fig. 14. The plot shows, on log-log scale, the variation of T with lattice size L, for given 
structure parameter 5 = 0.0001. 

Fig. 15. Number n(s) of structures is plotted, on log-log scale, against their size s, for two 
dimensional lattice with size L x L = 100 x 100, for different W. Parameters taken 
are e = 0.6, /i — 4, and 5 = 0.0001. We have used open-boundary conditions. 3000 
transients are discarded. Data are obtained for 20000 iterates per initial conditions 
and with 4 initial conditions. 

Fig. 16. The plot shows variation of average size s of the two dimensional structures with noise 
W, for the parameters as stated in Fig. 15. Error bars indicate the standard deviations 
at the data points. 

Fig. 17. Distribution of n(r) is plotted for two dimensional CML, on log-linear scale, against 
the lifetime r of the structures, for system parameters as in Fig. 15. 

Fig. 18. Plot shows variation of r with noise W , for two dimensional CML. 

Fig. 19. Distribution of the number n(S 2 ) of STS is plotted, on log-log scale, against the size 
S 2 , for different noise strength values W, for a two dimensional CML, parameters being 
same as in Fig. 15. 

Fig. 20. Variation of average size S 2 is shown against noise strength value W for a two dimen- 
sional CML. 
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